Bubble nucleation in stout beers 
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Bubble nucleation in weakly supersaturated solutions of carbon dioxide — such as champagne, sparkling 
wines and carbonated beers — is well understood. Bubbles grow and detach from nucleation sites: gas pock- 
ets trapped within hollow cellulose fibres. This mechanism appears not to be active in stout beers that are 
supersaturated solutions of nitrogen and carbon dioxide. In their canned forms these beers require additional 
technology (widgets) to release the bubbles which will form the head of the beer. We extend the mathemat- 
ical model of bubble nucleation in carbonated liquids to the case of two gasses and show that this nucleation 
mechanism is active in stout beers, though substantially slower than in carbonated beers and confirm this by 
observation. A rough calculation suggests that despite the slowness of the process, applying a coating of hollow 
porous fibres to the inside of a can or bottle could be a potential replacement for widgets. 

PACS numbers: 47.55.db, 64.60.qj, 82.60.Nh 



I. INTRODUCTION 

The production of bubbles in weakly supersaturated solu- 
tions of carbon dioxide is of great interest to the beverage in- 
dustry. Such solutions include many soft drinks and beers, as 
well as sparkling wines and champagne. While it has long 
been appreciated that spontaneous bubble formation in these 
liquids is strongly inhibited and thus that bubble formation 
can only occur at certain nucleation sites JH Hol| . it is only 
comparatively recently that the nature of these sites has been 
fully elucidated. In a series of papers, Liger-Belair and co- 
workers demonstrated that the most important nucleation sites 
are pockets of gas trapped in cellulose fibres Q (an example 
of type IV nucleation in the classification of Jones et al. yD) 
and developed a mathematical model of the growth and de- 
tachment of these bubbles @|, (a complementary model mak- 
ing slighly different assumptions was developed by Uzel et 
al. @]). 

While most beers are carbonated, there are advantages to 
using a mixture of nitrogen and carbon dioxide in beers, as 
is done in a number of stouts. (Hereafter, the term 'stout' 
will be used to refer to a beer containing a mixture of dis- 
solved nitrogen and carbon dioxide.) These advantages in- 
clude lower acidity in the beer leading to an improved taste; 
and smaller bubbles giving a creamy mouthfeel and a long 
lasting head dlH]. These beers are interesting scientifically 
because they show interesting fluid dynamical phenomena 
such as roll waves H] and sinking bubbles iflTll . Also of sci- 
entific interest is the technology used to create the head in the 
canned products. 

Pouring a carbonated beer from the can into a glass is 
enough to generate the head. This is not the case for stouts. 
Foaming in canned stouts is promoted by a widget: a hollow 
ball containing pressurised gas. When the can is opened, the 
widget depressurises by releasing a gas jet into the beer. The 
jet breaks up into tiny bubbles which are carried throughout 
the liquid by the turbulent flow generated by the gas jet and by 
pouring the beer into a glass. Dissolved gasses diffuse from 



TABLE I. Values of parameters used in this work. 



Parameter 


Value 




Reference 


r 


6.00 x 10- 6 m 




16] 


A 


14.00 x lO^m 




16] 


7 


47.00 x lO^Nm" 1 




HI 


Dt 


1.40 x 10- 9 m 2 s- 1 






D 2 


2.00 x 10- 9 m 2 s- 1 






#i 


3.4 x 10~ 4 molrrr 


i N -i 




H 2 


6.1 x 10~ 6 molirT 


i N -i 




T 


293 K 






Po 


1.00 x 10 5 Pa 






Pi 


0.80 x 10 5 Pa 




17] 


Pi 


3.00 x 10 5 Pa 
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the liquid into the bubbles which rise to the surface of the beer 
to form the head. 

In this paper we extend the mathematical model of bubble 
formation in carbonated liquids developed by Liger-Belair et 
al. 16[] to the case of two dissolved gasses and use it to inves- 
tigate two questions: 

• Why do stout beers require widgets? Is the bubbling 
mechanism described by Liger-Belair et al. completely 
inactive in stout beers or merely very slow? 

• Could an alternative to the widget be developed by coat- 
ing part of the inside of the can by hollow fibres? 



II. MATHEMATICAL MODEL 

In this section we develop a mathematical model of the rate 
of growth of a gas pocket in a cellulose fibre for the case 
in which there are two dissolved gasses: nitrogen and car- 
bon dioxide. Once a gas pocket reaches a critical size (when 
it reaches an opening of the fibre) it rapidly forms a bub- 
ble outside the fibre, leaving behind the original gas pocket. 
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FIG. 1. Geometry of a gas pocket trapped in a cellulose fibre. 



Since bubble formation and detachment is much faster than 
the growth of the gas pocket, the rate at which bubbles are 
nucleated can be deduced from the rate of growth of the gas 
pocket H]. 

The geometry of a gas pocket in a cellulose fibre is shown 
in Fig. [TJ Dissolved gasses in the fluid diffuse into the bub- 
ble through the walls of the cellulose fibre and through the 
spherical caps at the ends of the gas pocket. The rate at which 
this occurs is determined by the surface area, the diffusion 
constant and a diffusion length scale. The diffusion constants 
used to calculate the fluxes of carbon dioxide and nitrogen 
through the spherical caps are the diffusion constants in free 
fluid: Di and D 2 . The relevant diffusion constants for flow 
through the cellulose walls are D\± and D 2 ±. NMR mea- 
surements show that for carbon dioxide D\± ps 0.2£>i 10]. 
We assume the same relationship holds between D 2 ± and D 2 - 
The diffusion lengthscale A was measured experimentally for 
carbon dioxide [,6], again we assume that this value is also 
valid for nitrogen diffusion. 



tion: 
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Ac 2 =H 2 [P 2 
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where Pi is the partial pressure of dissolved carbon dioxide, 
P-2 is partial pressure of dissolved nitrogen, P B is the pressure 
in the gas pocket given by the Laplace law, Pq is atmospheric 
pressure and 7 surface tension. 

These equations can be non-dimensionalised by using the 
scales 

iV sca ie = l^ 2 Pb ^ 3 « 3.22 x lCT 13 mol, (7) 



D 2± RT 



t scale — 



rP B \ 



2D 2± H 2 P 2 RT 



2.73 s, 



(8) 



to introduce dimensionless variables Nx, N 2 and t (by as- 
sumption D 2 /D 2 ± = Di/Di±). The dimensionless equa- 
tions are 



aiiVi 



e ^ = (1 + JVl+W2 ,(i Ni + Ni 

dN 2 , ( a 2 N 2 



(9) 
(10) 



Using values typical of stouts, given in Table[I] the dimension- 
less parameters are 



D 2 H 2 P 2 

£ = ibtA w °- 096 ' 

a x = ^ w 1.45, 
"1 

Pb 

a 2 = -f w 0.39. 
"2 



(11) 
(12) 

(13) 



In this model the rate of change of the numbers of carbon 
dioxide (N*) and nitrogen (N 2 ) molecules in the gas pocket 
are given by 

^=4^Dx^+2,rzDx^, (1) 
at* A A 

^=^+^, (2) 
at* A A 

where asterisks indicate dimensional variables that will be 
non-dimensionalised later. 

Using Henry's law, Laplace's law and the ideal gas equa- 



III. ASYMPTOTIC SOLUTION 



Equations [9] and [10] cannot be solved directly. They can, 
however, be solved in two asymptotic limits: e <C 1 and 
Ni + N 2 1. The former limit does not produce particu- 
larly accurate results but the analysis of this limit helps us to 
interpret the results from taking the second asymptotic limit. 
The results from taking the second asymptotic limit are more 
accurate but harder to understand intuitively. 
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TABLE II. Numerical values of the parameters in equations |2T1 
andEU 



L £11 tl 1 1 1L LCI 




an 




a\i 


0.836 


a 2 \ 


-0.145 


022 


0.548 


n 


0.161 




0.468 


^scale 


0.439 s 




1.278 s 



A. First asymptotic limit: e <C 1 



Taking the limit in which the small parameter e « 0.1 is 
zero, equation[9]becomes an algebraic equation 



= 1 



aiNx 



N 1 +N 2 

which can be substituted into equationfTol 



dN2 
dt 



Oil + Oi 2 — OL\OL2 



OL\ — 1 

This equation is solved by 



No 



Oil + a 2 ~ OL\(X2 

ai 



. t \ oti + a 2 - a\a 2 

N 2 = ,4 exp - . 

t I a\ 



(14) 



(15) 



(16) 



where A is a constant of integration and r a dimensionless 
time constant describing the timescale of growth of the gas 
pocket in this approximation 



ai — 1 



Oi\ + OL 2 — OL\OL 2 

Ttscak = 0.954 s. 



0.35, 



(17) 
(18) 



Physically this approximation corresponds to assuming that 
diffusion of carbon dioxide is infinitely fast, and thus the par- 
tial pressure of carbon dioxide in the gas pocket is always 
equal to the partial pressure of carbon dioxide in solution. Ob- 
viously this approximation is only valid if the partial pressure 
of carbon dioxide is less than the gas pocket pressure, oth- 
erwise equation [141 has no physical solutions. Note that this 
approximation will underestimate r since it assumes carbon 
dioxide diffusion is infinitely fast. 



B. Second asymptotic limit: JVi + N 2 ^> 1 



In the limit Ni + N 2 3> 1 equations l9l andflOlbecome 



dNi 

~~dT 

dN 2 

— = N 1 + (l-a 2 )N 2 . 



.^1ni + ^, 

e e 



These equations have two independent solutions 



N\ = Aan exp — 



and 



Ban exp 



t 

Tl 
t 

T 2 



N 2 = Aa 2 i exp - 



t 

Tl 



N 2 = Ba 22 exp 



(19) 
(20) 

(21) 

(22) 



where A and B can be chosen independently to satisfy initial 
conditions. The numerical values of the other parameters are 
given in Table ITT1 

The analysis of the e < 1 case allows us to interpret these 
two solutions. The first solution, equation |2U decays ex- 
ponentially with a small timescale. This corresponds to the 
rapid establishment of the (dynamic) equilibrium concentra- 
tions (or partial pressures) of CO2 and N2 within the gas 
pocket (assumed instantaneous in the previous analysis). The 
second solution, equation [22] shows exponential growth with 
a longer timescale. This describes the steady state growth of 
the gas pocket at a fixed concentration ratio of CO2 to N2. 
The timescale of this process describes the timescale of bubble 
production. This analysis produces a longer estimate of that 
timescale than the previous analysis. This is because, previ- 
ously, diffusion of CO2 was assumed to be instantaneous. As 
the numerical results described in the next section show, the 
e <C 1 limit underestimates the correct timescale, while the 
Ni + N 2 ^> 1 analysis gives a good estimate. 



IV. NUMERICAL SOLUTION 

A full solution of the dimensionless equations can be ob- 
tained by numerical integration. A fourth order Runge-Kutta 
scheme with a timestep of 10~ 3 was used to solve equations|9] 
and [10] with initial conditions Ni = 0, N2 = 0.5. The differ- 
ential equations were solved over the interval < t < 10. 
The result for N = Ni + N 2 for 5 < t < 10 were fit- 
ted to an exponential curve giving a dimensionless bubble 
growth timescale of r = 0.47 corresponding to a dimensional 
timescale of rt sca i e = 1.28 s, in agreement with that predicted 
from the analysis of the Ni + N 2 ^> 1 case. This can be com- 
pared with the value for carbonated liquids at the same total 
pressure: 0.079 s. Fig. [2] shows the results of the numerical 
simulations over the interval < t < 1, 

In conclusion, these analytic and numerical results suggest 
that the mechanism of bubble formation described by Liger- 
Belair et al. is potentially active in stout beers but acts much 
more slowly than in carbonated drinks. 



4 



ln(iVi+iV 2 ) 
5 T 



(a) 




1 t 



N 1 +N 2 


0.8 -r 




0.6 - 




0.4 - 


f 


0.2 - 







h 



+ 



4- 



+ 



H t 



0.2 0.4 0.6 0.8 1.0 



FIG. 2. Results of the numerical solution of equations [9] and [10] 
The black line shows the numerical solution, the grey line shows the 
N\ + N2 S> 1 limit, and the dashed black line shows the e -C 1 
limit, (a) Rate of growth of the gas pocket, (b) Evolution of the 
concentration of CO2 in the gas pocket. 



V. EXPERIMENTAL CONFIRMATION 

To confirm experimentally that cellulose fibres could nu- 
cleate bubbles in stout beer, we observed a canned draught 
stout in contact with cellulose fibres (taken from a coffee fil- 
ter). Before opening the can, we made a small hole in the 
can to slowly degass the widget. This prevented it from foam- 
ing, which would have removed the dissolved gasses from so- 
lution. Using a microscope we observed that bubbles were 
indeed nucleated from the cellulose fibres but at a relatively 
slow rate. Fig.[3]shows bubbles nucleated by a gas pocket: the 
three parts of Fig.[3]are frames taken from a movie. 

Fig.|4]shows the growth of the gas pocket shown in Fig. [3] 
The figure was constructed from the same movie of the bub- 
bling process used for Fig. [3] Two hundred frames, corre- 
sponding to 8 s were extracted from the movie and rotated so 



that the fibre shown in Fig. [3] was vertical. From each frame 
the same column of pixels, passing through the centre of the 
fibre, was extracted and those columns placed side by side to 
construct a new figure: Fig. [4] This figure shows the evolution 
of the gas pocket: its slow growth (as predicted by the model) 
and then its rapid loss of gas to form an external bubble (as 
assumed by the model). 



VI. WIDGET ALTERNATIVES 

The model developed above allows us to investigate the fea- 
sibility of an alternative foaming strategy for stout beers in 
cans and bottles in which a coating of hydrophobic fibres on 
the inside of the can is used to promote foaming. A typical 
pouring time for a stout beer is 30 s. In this time about 10 8 
postcritical nuclei must be released. A single fibre produces 
one bubble every 1.28 s. Therefore about 4.3 x 10 6 fibres are 
needed. If each fibre occupies a surface of area A 2 then the 
total area that must be occupied by fibres is 8.3 x 10~ 4 m 2 
equivalent to a square with edge length 2.9 cm. This indicates 
that such an approach may be possible. 



VII. CONCLUSIONS 

A model of bubble formation in carbonated liquids has been 
extended to the case of liquids containing both dissolved ni- 
trogen and carbon dioxide. Taking values typical of stout 
beers shows that bubble formation by this mechanism does 
occur but at a substantially slower rate. This is consistent with 
the observation that widgets are needed to promote foaming 
in canned stouts. The possibility of replacing widgets with 
an array of hollow fibre nucleation sites was investigated and 
shown to be potentially feasible. 
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FIG. 3. Bubble nucleation in stout from a cellulose fibre taken from a coffee filter. The scalebar is 50 urn in each figure. The arrows show: 
1 a gas pocket in the fibre which nucleates a bubble; 2 the bubble fed by gas from pocket 1 ; 3 a bubble nucleated on the outside of the fibre, 
(a) The air pocket (1) has reached maximum size, (b) The air pocket (1) has created bubble (2). (c) Bubble (2) has visibly detached from the 
fibre, air pocket (1) is starting to refill with gas. (a), (b) and (c) are frames from a movie, (b) is 80 ms after (a), (c) is 520 ms after (a). 
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FIG. 4. Growth of the gas pocket. This figure shows the growth 
of the gas pocket within the cellulose fibre shown in Fig. [5] As in 
Fig.|3]dark colours correspond to gas and light colours to liquid. The 
columns of pixels from the frames corresponding to parts (a), (b) and 
(c) of Fig. [3] are indicated. The location of the gas pocket indicated 
by a (1) in Fig.[3]is shown. 



